Influence of antisymmetric exchange interaction on quantum tunneling of 
magnetization in a dimeric molecular magnet Mug 
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We present magnetization measurements on the single molecule magnet Mng revealing various 
tunnel transitions inconsistent with a giant spin description. We propose a dimeric model of the 
molecule with two coupled spins S — 6 , which involves crystal field anisotropy, symmetric Heisen- 
berg exchange interaction and antisymmetric Dzyaloshinskii-Moriya exchange interaction. We show 
that this simplified model of the molecule explains the experimentally observed tunnel transitions 
and that the antisymmetric exchange interaction between the spins gives rise to tunneling processes 
between spin states belonging to different spin multiplets. 

PACS numbers: 75.45. +j, 75.50.Xx, 75.60.Ej, 75.60.Jk, 75.75.+a 



Single-Molecule Magnets (SMMs) have been studied 
intensively in recent years because of the unique crossover 
between classical an quantum physics [H El [3j H] . The 
macroscopic observation of quantum phenomena such as 
tunneling between different spin states or quantum in- 
terference between tunneling paths give the possibility 
of studying in detail the quantum mechanical laws in 
nanoscale molecular systems and also might provide sub- 
stantial information concerning the implementation of 
spin based solid state qubits |51 IHl 13 H] • 

During the last ten years the spin system of SMMs has 
mainly been described by a single macroscopic spin and 
the associated tunneling processes were transitions inside 
a multiplet with total spin S, i.e. transitions that con- 
serve the total spin S of the molecule [5] [TUJ [H] . Recent 
developments in the field of molecular magnetism go be- 
yond this giant-spin approximation [131 HH [HI [IB] . 
When describing the molecule as an object composed of 
several superexchange-coupled spins Si the total spin S 
of the molecule is not fixed, but several multiplets with 
different total spin S appear and as a consequence the al- 
lowed tunnel transitions and relaxation paths of the spin 
system increase considerably. The associated tunnel pro- 
cesses between different spin states in this multi-spin de- 
scription do not need to conserve the total spin S of the 
molecule. Recently Carretta et al. showed evidence of 
this quantum superposition of spin states with different 
total spin length in the molecule magnet CryNi by inelas- 
tic neutron scattering (INS) [T3j. In fact, when introduc- 
ing an antisymmetric exchange coupling (Dzyaloshinskii- 
Moriya interaction) between the spins Si that compose 
the molecule, the superposition of a symmetric and an 
antisymmetric spin state becomes possible. The associ- 
ated tunneling process and quantum interference effects 
of different tunneling paths have been observed recently 
in a Mni2-based molecular wheel 

In this Brief Report we report the observation of quan- 



tum tunneling between spin states with different total 
spin S in a Mn-based SMM having presently the high- 
est anisotropy barrier of 89 K [17] . A theoretical model 
is proposed that describes the molecule as an exchange- 
coupled system of two separated spins Si — 6. The 
experimentally obtained tunnel splittings that uses the 
Landau-Zener method of various symmetric and antisym- 
metric tunnel transitions are compared to this theoretical 
model. 

The SMM has the chemical formula [Mnl"02{Et ~ 
sao)ei02C'Ph{Me)2)2{EtOH)e] and will be called briefly 
Mtiq ^TT}. The six Mn atoms, each having a spin Si = 
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FIG. 1: (color online) Zeeman diagram of the dimeric 
molecule Mng using the longitudinal anisotropy constant 
D — 1.28 K, an isotropic Heisenberg exchange interaction 
J = 0.8 K and g — 1.99. Due to the exchange interaction 
some excited spin multiplets are located only a few Kelvins 
above the ground state. The inset shows a simplified model 
of the magnetic core of the Mne molecule. Two ferromagnet- 
ically coupled spin triangles, each having a total spin S = 6, 
form the dimeric molecule. 
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2, form the core of the molecule and they are strongly 
superexchange coupled and act as a macroscopic spin S — 
12 at low temperature. Recent work by Carretta et al. 
employing INS shows evidence of very low lying excited 
spin multiplets in Mng resulting in the breakdown of the 
giant spin model [TB]. In contrast to Carretta et al. we 
propose to describe the molecule by two superexchange 
coupled spin triangles, each of them being described by a 
rigid total spin S = 6 (see inset of Fig. 1). This molecular 
dimer description is in very good agreement with the INS 
measurements and simulations shown by Carretta et al. 
This simplified model reproduces very well the low-lying 
spin multiplets and gives the advantage of a small Hilbert 
and parameter space compared to the description of the 
Mug molecule in Ref. [IB] . 

Each of the two ferromagnetically coupled spins of the 
molecular dimer S'l = 5*2 = 6 can be described by the 
spin Hamiltonian : 

- -DiS^f + 0(4) - .9//BMo5. • H (1) 

where Sf, Sf and are the vector components of the zth 
spin operator, g = 1.99 is the gyromagnetic factor and 
/iB is the Bohr magneton. [T7] The first term describes 
the uniaxial anisotropy of the molecule with longitudinal 
anisotropy parameter D and the second term contains 
fourth order crystal field anisotropy terms. The last term 
is the Zeeman interaction of the spin 5*^ with an external 
magnetic field H. 

The exchange interaction of the two halves of the 
molecule can be described by 

Hex = -J Si ■ S2 + Di2 ■ (Si X ^2) (2) 

where the first term describes the isotropic Heisenberg 
exchange interaction with exchange constant J and the 
second term is an antisymmetric Dzyaloshinskii-Moriya 
interaction between the two spins. 

Exact diagonalization of the total spin Hamiltonian 
Ti. = Til + 7^2 + T^cx leads to the energy spectrum shown 
in figure 1. The lowest lying spin states belong to the 
S" = 12 multiplet. Due to the ferromagnetic exchange 
the first excited spin multiplet jS* = 11, Ms = ±11) is 
located at about 25 K above the ground state doublet 
\S — 12, Mg — ±12) in zero magnetic field. 

In the following we will discuss the different level cross- 
ings not in the eigenbasis of the total spin of the molecule 
\S,Ms). As the total spin of the molecule may fluctuate 
we chose the more convenient eigenbasis of the two sin- 
gle spins of the molecule \Si,mi) (g) 15*2,771.2) = |™i,m2). 
The groundstate doublet can be expressed as |12, ±12) = 
I ±6, ±6) and the first doublet of the first excited multi- 
plet reads 111, ±11) = ^(|±6,±5)-|±5,±6)). The low- 
est lying spin eigenstates belonging to the ground state 
multiplet S = 12 are symmetric in respect to a permu- 
tation of the two spins in the product base |toi) (8) \m2), 
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FIG. 2: (color online) Magnetization measurements for dif- 
ferent field sweep rates and two transverse magnetic fields 
jj,oHx = 3.4 T (solid lines) and I-LqHx = 4.1 T (open symbols). 
The sample was first saturated in a large negative magnetic 
field and then ramped at constant sweep rate to positive field. 
All measurements were done at low temperature T = 100 mK. 



whereas the eigenstates of the 5 = 11 multiplet are an- 
tisymmetric. When we look at the probability to tunnel 
from one spin state to another, we see immediately that 
most of the terms in the Hamiltonian Ti. are symmetric 
and therefore only provide coupling between symmetric 
spin states. The only antisymmetric term in the Hamilto- 
nian is the Dzyaloshinskii-Moriya exchange interaction, 
and as a consequence this term can provide a coupling 
between a symmetric and an antisymmetric spin state, 
i.e. this term couples spin states of the ground state 
multiplet S* = 12 and the first excited multiplet S* = 11. 

The magnetic measurements on a single, /im-sized 
crystal were carried out in a dilution refrigerator employ- 
ing a vector magnet system and a Hall sensor [TF| . The 
easy axis of the magnetization of the crystal was paral- 
lel to the z-direction of the applied magnetic field. The 
tunnel splittings of the anticrossings were determined fol- 
lowing the Landau-Zener technique [51 [TOl EUl 1^ . 

Figure 2 shows some magnetization measurements at 
low temperature for large transverse magnetic fields at 
different field sweep rates. Equally spaced and very pro- 
nounced steps of the magnetization appear at approxi- 
mately ^j.qH^ « 0, 0.45, 0.9, 1.35 and 1.8 T. In between 
these tunnel transitions, we observe a fine structure of 
smaller steps, which occur at approximately ^qH^ ~ 1.2 
and 1.65 T. 

Figure 3 shows the derivatives of the magnetization 
curves of figure 2 as well as the corresponding Zeeman 
diagram with the lowest energy levels. The main steps of 
magnetization, equally spaced by AfioHz « 0.45 T, can 
be explained in the framework of a giant spin approxi- 
mation, when describing the molecule by a collective spin 
S = 12. We checked that the fine structure is not due to 



3 




10"' 



0.5 1.0 1.5 

Magnetic field (T) 

FIG. 3: (color online) (a) Zeeman diagram with tunnel tran- 
sitions between symmetric spin states within the ground state 
multiplet and tunnel transitions from the ground state mul- 
tiplet to spin states belonging to excited spin multiplets. (b) 
The derivative plot of the magnetization curves of figure 2 
shows various peaks due to tunnel transitions. In between 
the main tunnel transitions between states belonging to the 
ground state multiplet we observe a fine structure of addi- 
tional tunnel transitions involving excited spin multiplets. 



spin-spin cross-relaxation [2U]. However, the fine struc- 
ture in the magnetization steps is related to excited spin 
multiplets. These steps can be understood when consid- 
ering excited spin multiplets in the multi-spin approach. 
The tunnel transition at hqHz ~ 0.45 T involves the 
symmetric eigenstates |— 6,— 6) and ^(|6, 5) -I- |5, 6)). In 
fact, in between the main, equally spaced tunnel tran- 
sitions several avoided level crossings appear involving 
excited spin multiplets (as shown by the blue dots in 
figure 3a). As an example, the avoided level cross- 
ing at fJ-oHz ~ 0.75 T involves the symmetric eigenstate 
|— 6,— 6) and the antisymmetric eigenstate ^(|6, 5) — 
|5,6)). The tunnel process at /.toiJ^ ~ 1.2 T involves 
the symmetric eigenstate |— 6,— 6) and the antisymmet- 
ric eigenstate ^(|6,4) — |4, 6)). The observed avoided 
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FIG. 4: (color online) Tunnel splittings A as a function of the 
transverse magnetic field for different level anticrossings. The 
longitudinal magnetic field hqH^ was swept over an avoided 
level crossing at a constant sweep rate dH^/dt — 68 mT/s and 
with fixed transverse magnetic field. The tunnel splittings 
were obtained by applying the Landau-Zener formula. The 
inset shows the time dependence of the magnetization in the 
saturated sample when sweeping several times over the level 
anticrossing. 



level crossings in our experiments allow us to determine 
the longitudinal anisotropy parameter D — 1.28 K and 
the isotropic exchange constant J — 0.8 K. 

Figure 4 shows the tunnel splittings A of the differ- 
ent level anticrossings within the ground state multiplet 
(at fioHz ~ 0.45, 0.9 and 1.35 T) and the ones involving 
excited spin multiplets (at /io-ffz ~ 1.2 and 1.65 T) as 
a function of the transverse magnetic field fioHx- Note, 
that the tunnel splittings of the two antisymmetric level 
anticrossings around figHz ~ 1.6 T could not be stud- 
ied seperately as they are too close. In order to deter- 
mine the tunnel splittings the longitudinal magnetic field 
was swept over a level anticrossing with fixed sweep rate 

dt 



— 68^ and fixed transverse magnetic field noH^ 
and the probability of tunneling from one state to the 
other was measured by means of the magnetization de- 
crease of the saturated sample. The tunnel probability 
Pm.m' between two spin states m and m' is given by the 
Landau-Zener formula 



exp 



fiodHjdt 



which allows calculation of the tunnel splitting of the 
avoided level crossing Aj„ „i' when Pm,'m' <C 1 UHl [201 
IHJ HI]. Note, that P™,„/ < 1 is not fulfilled for very 
high transverse magnetic fields and therefore the experi- 
mentally obtained tunnel splittings are only estimates of 
a lower bound of A^.m'- The experimentally obtained 
tunnel splittings lie in the range of 10^^ K for all the 
observed transitions and they rapidly increase when ap- 
plying a transverse magnetic fields fioHx > 3 T. 
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We found that the tunnel splittings of the anticross- 
ings between symmetric states are mainly determined 
by the symmetric spin operators such as the second 
and fourth order anisotropy terms or the Heisenberg ex- 
change interaction. However the splitting between a sym- 
metric and an antisymmetric spin state is given by the 
matrix element involving antisymmetric spin operators, 
i.e. in the framework of our model the antisymmet- 
ric Dzyaloshinskii-Moriya interaction. When analyzing 
the tunnel splitting between a symmetric and an an- 
tisymmetric spin state we can get an estimate of the 
Dzyaloshinskii-Moriya interaction parameter Di2. Fur- 
ther on, the magnitude of the tunnel splitting between 
symmetric spin states can be used to fix the parameters 
D, J and possible fourth order parameters. 

Numerical simulations of the tunnel splittings by exact 
diagonalization of the above defined Hamiltonian show, 
that the isotropic exchange and weak higher order trans- 
verse anisotropy terms together with a transverse mag- 
netic field comparable to the one used in the experiments 
give rise to tunnel splittings between states of the ground- 
state multiplet on the order of 10""^ K. The magnitude of 
the tunnel splittings is well reproduced when introduc- 
ing a weak fourth order spin operator term as proposed 
by Carretta et al. |16j The strong increase in the tun- 
nel splitting Am,m' for large transverse magnetic fields is 
also well reproduced in the framework of this model. 

The large tunnel splittings between symmetric and an- 
tisymmetric spin states cannot be reproduced by any 
symmetric spin operator such as second and fourth order 
crystal field anisotropy terms or the Heisenberg exchange 
interaction. However, the antisymmetric Dzyaloshinskii- 
Moriya exchange interaction can provide the quite large 
coupling between the symmetric and antisymmetric spin 
states. Numerical simulations show that a vector of 
the Dzyaloshinskii-Moriya interaction D12 with compo- 
nents Dx = Dy = Dz = 10 mK gives rise to tunnel 
splittings of the order of the experimentally observed 
ones. In particular, the calculated tunnel splitting at 
HqHz ~ 0.75 T turns out to be at least one order of mag- 
nitude smaller than the one at fJ-oH^ ~ 1.2 T. This is con- 
sistent with our experiments, as we did not observe any 
clear and pronounced step in the magnetization curves 
at fioHz ~ 0.75 T. The corresponding tunnel splitting is 
theoreticaly - with the parameters given above and trans- 
verse magnetic fields below 4 T - smaller than 10~* K and 
therefore too small to be measured with our experimental 
technique. 

In conclusion, we presented magnetization measure- 
ments on the SMM Mug revealing various tunnel tran- 
sitions, which are forbidden in the framework of a giant 
spin approximation. We propose to describe the Mng 
SMM as a molecular dimer of two coupled spins S = 6. 
The introduction of an antisymmetric exchange interac- 
tion leads to the superposition of spin states with dif- 
ferent spin length. This superposition of spin states be- 



longing to different multiplets leads to additional tunnel 
transitions which are observed in our experiments and 
are in perfect agreement to our theoretical model. This 
multi-spin description goes far beyond the standard gi- 
ant spin approximation and is capable of explaining the 
experimentally observed tunnel transitions. This dimeric 
model of the molecule is confirmed by numerical calcu- 
lations of the positions and the magnitude of the tunnel 
splittings which are consistent with the experimental re- 
sults. 
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